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Abstract 

Based on empirical financial time-series, we show that the "silence-breaking" 
probability follows a super-universal power law: the probability of observing 
a large movement is inversely proportional to the length of the ongoing low- 
variability period. Such a scaling law has been previously predicted theoreti- 
cally [ij], assuming that the length-distribution of the low- variability periods 
follows a multiscaling power law. 
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1. Introduction 

The power laws and scaling behaviour are present in numerous aspects of 
human societies and in the nature. One of the first promoters of the concept 
of the power law was Vilfredo Pareto (cf. 0|), who studied the wealth dis- 
tribution in different societies. Further, Harvard linguistics professor George 
Zipf, again a representative of social sciences, observed that only few words in 
English language are used very often, and most of the words are used rarely, 
(cf. [^). Nowadays, physicists are very used to the power-laws, which, how- 
ever, are sometimes somewhat counter-intuitive. Indeed, the presence of a 
power-law means that there are some representatives of a population, which 
are very different from the typical members of that population. For exam- 
ple, as the result of a social evolution, the wealth of a single individual can 
qualitatively change the wealth of a large community. This is completely 
different from the biological evolution: e.g. the weight of a single living 
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creature makes always only a tiny contribution to the net biomass of the 
corresponding (non-small) population. 

The presence of a wide spectrum of power laws in finances can be ascribed 
to the fact that socially interacting humans form large complex systems, 
which are characterized by self-organized criticality This makes finances 
(alongside with the turbulence) a fruitful polygon for studying various as- 
pects of scale-invariance. Indeed, various power-laws have been observed in 
financial time series since 1960-es, by B. Mandelbrot (cf. [11,01 and references 
therein). In 1999, Ausloos and Ivanova also reported the multifractality in 
financial time series (cf. 01). Around 1990-ies, the studies of scale-invariance 
in finances became more extensive, cf. 0, [l^, 11|, effectively creating a 



new branch of statistical physics — the econophysics. A recent overview of 
the progress in understanding the scaling and its universality in finances can 



be found in Ref. The aim of the current study is to contribute to the 
understanding of the origins of universality. Mathematically, our basic idea 
is very simple, nearly trivial. However, when dealing with the sources of 
universality, mathematically simple and robust models have better chances 
of describing reality, than complex and elaborate constructions; cf. the Oc- 
cam's razor. 

The attempt to successfully and systematically predict the direction of 
future movements of asset prices can be compared to the attempts of invent- 
ing the perpetum mobile. However, the attempt to characterize and predict 
the risk (or volatility) may offer significantly better results and therefore the 
volatility is one of the most-studied phenomena in Econophysics. 

One of the most challenging features of the volatility dynamics are the 
intermittently appearing extreme price movements, which are often accom- 
panied by overall increase of volatility over a certain time window. Tradi- 
tionally, such a behaviour has been described by multifractal spectra. The 
multifractal analysis is untoubtedly a powerful tool; however, due to the in- 
volved mathematical methods, it is not well-suited for practical applications 
of the prediction and optimisation of risks. This observation motivated us to 
introduce a complementary method of the length-distribution analysis of the 



low-variability periods [13|, llJ, Il5||. The low- variability periods are defined 
as consequent time periods, during which the price changes of the observed 
asset (as compared to the local average over a sliding window of width w) 
remains under a pre-set threshold 6. The illustration of the method is shown 
in Fig. [H This method (with certain modifications) was developed indepen- 
dently within different contexts, and put under extensive tests by several 
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research groups. So, the effects of the long-term memory and clustering of 



extreme events in various time series (cf. [l6|, ll7|, ll8|, ll9|, l26||), are, in fact, 
closely related to (and in a certain sense covered by) the low-variability pe- 
riod analysis. The same applies to the studies of the time intervals r between 



volatilities which are above a threshold q [20|, l2l|, l22|, l23|, \24, l25l|). 



Our previous studies 15 1 and [l| have shown that 



• financial time-series are typically characterized by multi-scaling be- 
haviour of the low-variability periods (power law can be observed for 
a certain range of the parameters 6 and w, and the scaling exponent 
depends on these parameters); 

• theoretically, a multi-fractal time-series follows also a multi-scaling be- 
haviour of the low- variability periods; the scaling exponent can be ex- 
pressed via the multifractal exponents; 

• as compared to the multifractal analysis, the analysis of the low- variability 
periods is easier to implement, has higher resolution of time-scales, and 
the results of the analysis can be interpreted more straightforwardly. 

These findings agree well with the above cited independent studies. 

The market fluctuations have been also modelled as Levy flights and con- 
tinuum time random walks (CTRW); therefore it is also of interest to mention 
that in the case of uncorrelated Levy flights, there is no power-law for the 
length-distribution of the low-variability periods: length-distribution decays 
exponentially; the same applies to Gaussian time-series, e.g. (non)persistent 
random walks. Meanwhile, in the case of CTRW, there is a mono-scaling 
behaviour of the low-variability periods: the scaling exponent does not de- 
pend on the parameters 6 and w, and is defined by the exponent of the 
waiting-time distribution (27| . 

In addition to the above listed results, we have shown [H that the very 
presence of a power-law for the probability distribution function of the low- 
variability segment lengths (even if observed for a narrow range of the param- 
eters 6 and w) bears interesting consequence: the probability p of observing 
a large movement (exceeding the threshold parameter 6) in the time series 
during the next period of duration w, is inversely proportional to the length 
of the on-going low-variability period (i.e. to the time elapsed since the most 
recent large movement; note that by definition, p < 1, and p = 1 would 
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Figure 1: Variability A of the DAX index and the respective low-variabiHty periods for 4 
years (2003-2007), using w — 1 day (a) and w = 10 days (b). For a pre-fixed threshold 
level S — 2a (where a is the standard deviation of the signal) , low- variability intervals of 
duration n are formed as the intervals corresponding to such graph segments, which lay 
entirely inside the gray area. The small graphs in top share the time axis with the bottom 
graph, and illustrate the fragmentation of them into low-variability intervals by plotting 
the interval index i versus time t. 
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imply that for the given parameters, the next price movement exceeds def- 
initely the threshold S). This super-universal scaling law has been derived 



independently by Bogachev, et al [18|, ll9|. Here we repeat the derivation 
briefly. 

First, we assume that the length of the periods is measured in the units 
of the window length w. Then, the probability of a large movement during 
the unit time is given by the ratio of (a) the number of those low- variability 
periods, the length of which is exactly n, Na = R{n) — R{n + 1), and (h) 
the number of those low-variability periods, which have length m > n, Nb = 
R{n). Here, R{n) denotes the cumulative length-distribution function of the 
low- variability periods with length equal to n [i.e. R{n) represents the number 
of low- variability periods with length greater or equal to n time units]. So, 
we can express the probability of "silence-breaking" (i.e. the probability that 
the next movement exceeding the threshold 5) as a function of the on-going 
"silent" period n: 

p{n) = [R{n)- R{n + l)]/R{n). (1) 

If n is large, the difference R{n) — R{n + 1) can be calculated approximately 
as — Upon applying the power-law R{n) = RQn~°'^^'^\ we arrive at 
Na ~ aRon~°'~^ and Nf, ^ R^n'". Bearing in mind that p{n) = Na/Nf,, the 
final result is written as 

p{n) ^ arT^ . (2) 

Since the scaling exponent of this law is independent of the parameters 6 and 
w (as well as of the scaling exponent a), it can be called super-universal. The 
equation ([2|) can be also denoted as the law of the silence-breaking probability, 
because it describes the probability of a large movement after a longer "silent" 
period (characterized by small movements). Recently, similar application of 
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the risk measurement was proposed by Wang, et al 

As demonstrated above, the mathematics behind this super-universal law 
is extremely simple. However, we do believe that the consequences of it are 
profound, and it allows to shed light into the origin of universality in the 
dynamics of market fluctuations. In fact, a very similar behaviour (with 
a similar mathematical origin) has been detected in the context of direct 
avalanches in self-organised critical systems (c.f. l28i]), and proven to be a 
useful tool in understanding the universality in burst dynamics. These argu- 
ments motivated us to put equation ([2]) under several tests, using both the 
data of real market fluctuations, as well as the surrogate data. The aim is to 
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clarify, how well is the super-universal law followed by the real-world time- 
series, for which the scaling laws are far from being perfect (and is for certain 
applications better described, e.g., by stretched exponentials and for 

which the finite-size effects can (possibly) mask the theoretical scaling be- 
haviour. In order to show the presence of the 1/n-law in relatively short 
time-series (with unavoidably poor statistics of long low- variability periods), 
an appropriate data analysis technique has been developed (see Section 12.11) 

2. Empirical tests 

Further we proceed to test empirically the probability of the "silence- 
breaking". As changes in financial markets are believed to obey power-laws 
((29|). we use several financial time series (cf. Table [I]) in order to test the 
equation Some of these data are presented visually in Fig. [2l We also 
test the equation ([21) with artificially generated time series (i.e. surrogate 
data with known properties, see Appendix), which obey a power-law. 



Table 1: The data used in empirical analysis (source of all the data series is Bloomberg) 



Abbr 


Description 


Calendar Period 


# of trading days 


SPX 


Standard & Poor's 500 Index 


30/12/27 - 31/10/08 


20305 


DAX 


The German Stock Index 


01/10/59 - 31/10/08 


12326 


NKY 


Nikkei 225 Stock Average 


05/01/70 - 31/10/08 


9588 


MXEA 


The MSCI Europe, Australasia and Far 
East Index 


31/12/71 - 31/10/08 


9610 


CAC 


CAC-40 Index of Paris Bourse 


09/07/87 - 31/10/08 


5374 


UKX 


FTSE 100 Index 


03/01/84 - 31/10/08 


6282 


MXWO 


MSCI World Index 


03/01/72 - 31/10/08 


9610 


INDU 


Dow Jones Industrial Average 


03/01/1900 - 31/10/2008 


27333 



We analysed the daily returns of eight different indices of the world stock 
exchanges (cf. Table [T]). Cumulative distributions R{n) for the S&P 500 
index for various thresholds 6 are shown in Fig. [3l the other indices follow 
a similar behaviour. The values of 6 were normalized to the standard devi- 
ation, which was calculated for the whole sample of every index. We found 
that R{n) was in a good accordance with the power-law for a broad range of 
thresholds S (ranging from 3/4 to 4 standard deviations). For smaller thresh- 
olds, inertial range of the linear part in the log-log plot was too short for a 
meaningful scaling analysis. On the other hand, greater values of 6 produced 
relatively small total number of low-variability intervals, insufficient for the 
further analysis (cf. Table El). Therefore, the optimal range for the threshold 
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Figure 2: a) S&P 500 index; daily returns [measured in standard deviations] of S&P 500 
(b), German Stock Index DAX (c), and Nikkei 225 Index (d) 
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level was found to be between 1 and 2 standard deviations (cr) . Here we have 
adopted the value 6 = 2a. 
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Figure 3: Cumulative distributions of low-variability periods R{m > n) over the length n 
for the S&P 500 time series, using different threshold levels 

While the details of the functional relationship a{S, w) depend on the 
specific features of the (non-Gaussian) statistics of the returns, generic trends 
can be outlined as follows. Within the optimal range of the parameters, 
larger values of the threshold 6 lead to longer low- variability periods, and to 
a decrease of the overall number of the low-variability periods. The value 
of a is the (modulus of the) slope of the straight part of the R{n) curve in 
log-log graph (Fig. [3]), which typically decreases for such a process (i.e. with 
increasing S). Similar arguments lead to the conclusion that larger values 
of the sliding window width w lead to shorter low-variability intervals and 
hence — to larger values of a. 

The empirical tests of Eq. ^ were conducted as follows: (i) we count the 
number of low- variability periods R{n) with length equal to n; (ii) we calcu- 
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Table 2: Some characteristics of the distribution function R{n) for the S&P 500 index, 
using different values of the threshold parameter S and w = 1 day. N denotes the length 
of the longest low-variability period. 



5, in stddev 


R{1) 


i?,(10) 


N 


1/4 


4145 


2 


10 


1/2 


4394 


107 


59 


3/4 


3491 


294 


105 


1 


2523 


414 


215 


2 


713 


218 


949 


3 


274 


107 


1609 


4 


124 


58 


3086 


5 


67 


32 


6382 



late the probability p{n) = [R{n) — R{n + 1)]/ R{n) and {iii) plot it against 
n. In order to make the data presentation visually easier to understand, we 
modify step {iii) by replacing p{n) with p{n) x n. Thus, should the equation 
(l2l) hold, the plot of p{n) x n against n should yield a horizontal line. 

2.1. Finite-size effects 

Before proceeding to the results of the tests outlined above, it should be 
noted that the probability of silence-breaking cannot be expected to hold for 
all the values of n. We explain this claim in more details. 

1. Let us recall that in order to derive the scaling law (2), we approximated 
R{n) — R{n + 1) with — this is, however, legitimate only for large 
values of n. Therefore, one can expect that the super-universal power law 
(2) fails for small values of n. The surrogate data analysis confirms that 
this is, indeed, the case, see Fig. [H Using numerically generated very 
long data-series [which correspond to the distributions R{n) = Ron'""], 
we have calculated the quantity p{n)xn, and plotted it against n; in a full 
agreement with the theoretical expectations, it reaches asymptotically 
(for n 1) a plateau with a limit value p{n) x n a. 

2. The analysis of the empirical data shows that R{n) is not a smoothly 
decreasing function of n. Instead, for larger values of n, there are long 
plateaus with a constant value of R{n). Indeed, it is possible that for 
many values of n, there is no low- variability periods with the given length 
n. It is apparent that due to the finiteness of the data series, the higher 
values of n offer less statistics, i.e. the likelihood of observing a low- 
variability period of a length, equal to a fixed and very high value of n, 
is very low. Consequently, if we apply Eq. (1) directly to the empirical 
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Figure 4: The quantity p{n) x n is calculated for very long surrogate data series, and 
plotted against the length n; the results confirm the analytically predicted asymptotics 
p{n) X n = a for n » 1 

data series, the resulting empirical probability estimates p{n) follow a 
very singular behaviour: there are sharp peaks at those values of n, 
which happen to have a matching low-variability period; between these 
peaks, p{n) = 0. Evidently, such a simple-minded data analysis is not 
suited for testing the scaling law (2). 

Perhaps the most straightforward solution here would be to average 
the empirical probability estimates over a sliding window of a suitable 
length. However, this is certainly not the easiest and most elegant route, 
because it raises several technical issues (e.g. what is the best window 
width; it should variable and increase with n). Here, we have proceeded 
as follows. First, we approximate the step- wise experimental distribution 
law R{n) with a piece- wise linear function R{n). More specifically, for 
all these values of n (n = ni,n2, ■ ■ .n^) for which R{n + 1) 
we define R{n) — R{n); for all the other values of n we use a linear 
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interpolation: for n G (nj,nj+i), we define 

R{n) = R{ni){ni+i - n) + R{ni+i){n - ni)/{ni+i - rii). (3) 

Then, we can apply Eq. (1) to the smoothed distribution function R{n). 
3. Finally, even if we apply the smoothing procedure as described above, 
for very large values of n (of the order of the length N of the longest 
low- variability period), the statistical basis becomes so narrow that the 
function p{n) = [R{n) — R{n + l)]/_R(n) will fluctuate with a very large 
amplitude. Apparently, the statistical fluctuations around the theoret- 
ical power-law [Eq. ([2])] grow with n; these fluctuations set an effective 
upper limit for the empirical scaling range at Umax ~ N . 

Being guided by these considerations, we calculated the -R(n)-functions 
for all the empirical data series, and applied the linearization procedure ac- 
cording to Eq. dS]). The results for the S&P 500 are presented in Fig. [5l 
the graphs for our other empirical data series behave similarly. Note that in 
Fig. [HI we have limited the data to n < due to the reasons discussed 

in the previous pragraph. As predicted theoretically, the grahps fluctuate 
around a constant value [for the parameters of Fig. [5|, p{n)n ^ 0.8]. 

3. Discussions: the consequences of the power-law 

Using the financial time-series, as well as artificially generated (surro- 
gate) data series, we have shown that the "silence-breaking" probability fol- 
lows, indeed, a super-universal power law: the probability of observing a 
large movement is inversely proportional to the length of the ongoing low- 
variability period. The usefulness of this result is not limited to the financial 
risk analysis, and helps to understand the origins of universality in the scal- 
ing behaviour of complex systems; it can be also used to build a hierarchical 
risk prediction scheme (for a spectrum of events of different amplitude). 

On the one hand, our result "approves" the common human tendency to 
forget those extreme events, which have happened a long time ago. Indeed, 
we have shown that in any system producing a scale-invariant sequence of 
events of different amplitude, the probability of observing a large event on 
next day is inversely proportional to the time elapsed since the last large 
event. On the other hand, however, our finding does not void another uni- 
versal property of power laws: sooner or later, there will be an event of an 
even larger amplitude. 
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Figure 5: The quantity p{n) x n is plotted for the S&P 500 index {6 — 2a, a = 0.81) and 
for the numerically generated time series (of a similar length and with the same value of 
a) against the period length n. These results are in a good agreement with the theoretical 
expectation: for 1, the graphs fluctuate around the asymptotical value p{n) xn ^ a. 



So, in the case of stock markets, we can conclude that if we have not seen 
a ten-percent-or-larger movement for the last 10 000 days, the probability of 
such a movement on next day can be estimated as 10~^. However, one should 
not be driven to a delusion: due to the multi-scaling nature of the stock 
market movements, even a movement of 20% cannot be neglected; Nassim 
Nicholas Taleb (cf. [3o|) has compared such a delusion with the Russian 
Roulette with a revolver containing 500 bullet-holes: the fact that nothing 
has happened during the 100 first clicks on the trigger, does not make the 
game harmless: the lOP* click can be lethal. 
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A. Appendix 

In order to test the equation ([2]), artificial data with know statistical 
properties we generated as follows. To begin with, note that actual time- 
series are not needed; instead, it suffices to generate the sequence of the 
durations Tj {i = 1,2, . . .) of the low-variability periods. Indeed, our aim is 
to analyse the "silence-breaking" probability p(n) for a fixed set of parameters 
w and 6 (without intention to repeat the test for the same time-series with 
different parameter values). 

So, the total number of low- variability intervals was fixed to be L = 10^ for 
Fig. [Hand L = 10^ for Fig. [5l Further, for every index i G (0, L) the random 
power-law-distributed distributed number was obtained. The values of Tj 
were interpreted as the lengths of low-variability periods. Finally, all the L 
intervals were sorted according to their length and cumulative distribution 
function R{n) was calculated. 
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